Anti-GZK effect in UHECR spectrum 
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Abstract. In this paper we discuss the anti-GZK effect that arises in the framework of the diffusive propagation of Ultra High 
Energy (UHE) protons. This effect consists in a jump-Uke increase of the maximum distance from which UHE protons can 
reach the observer. The position of the jump is independent of the Intergalactic Magnetic Field (IMF) strength and depends 
only on the energy losses of protons, namely on the transition energy from adiabatic and pair-production energy losses. The 
Ultra High Energy Cosmic Rays (UHECR) spectrum presents a low-energy steepening approximately at this energy, which 
is very close to the position of the observed second knee. The dip, seen in the universal spectrum as a signature of the proton 
interaction with the Cosmic Microwave Background (CMB) radiation, is also present in the case of diffusive propagation in 
magnetic fields. 
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INTRODUCTION 

Recently a very interesting phenomenon, determined by 
the UHE proton propagation in IMF, has been found 
Jlll2l - It consists in a low energy steepening of the pro- 
ton spectrum that occurs at energy below 10^^ eV. This 
steepening is caused by an increase of the diffusive prop- 
agation time, that rapidly exceeds the age of the universe, 
and can be explained trough the diffusive propagation of 
UHECR in IMF. The position of the steepening energy 
Es — 1 X lO'** eV, is determined only by the proton en- 
ergy losses on CMB and coincides with a good accuracy 
to the position of the 2nd knee observed in the CR spec- 
trum 1 3]. In this paper we will discuss the main features 
of the diffusive propagation of UHECR in IMF, focus- 
ing our attention on the steepening energy scale E^. Be- 
fore entering the details of the diffusive propagation of 
UHECR protons let us review the main experimental ev- 
idences related to IMF. 

The presence of an IMF is still an open question, the 
most reliable observations of this field are based on the 
measurement of the Faraday rotation (RM) of polarized 
radio emission |4|]. The upper limit obtained with these 
measurements is RM < 5 rad/m^, it implies an upper 
limit on the IMF that depends on the assumed scale of 
coherence length. For instance, according to in the 
case of an inhomogeneous universe B/^ < 4 nG with a 
scale of coherence of about — 50 Mpc. In general, as 
follows from the observations of Faraday rotation, the 
magnetic field is high, of the order of 1 jiG with a co- 
herence length Ic = 1 Mpc, in clusters of galaxies and ra- 
dio lobes of radio galaxies |4|]. Apart from observations. 
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the IMF can be predicted, in principle, trough Magneto- 
hydrodynamics (MHD) simulations. The main ambigui- 
ties in these simulations are related to the assumed seed 
magnetic field and to the capability of simulations to re- 
construct the local Universe as we observe it (i.e. con- 
strained 1 6] and unconstrained simulations Q]). Unfortu- 
nately, because of these uncertainties, MHD simulations 
are not completely conclusive, there are at least two op- 
posite results in literature with predicted magnetic field 
in voids (filaments) that vary from 10-3 nG (10-' nG) 
161 up to 10-' nG(lOnG) |7]. 

While a direct evaluation of the IMF strength is still 
challenging, indirect informations about the UHECR 
propagation mode can be inferred from UHECR data. 
The analysis of the arrival directions of UHECR at ener- 
gies E > lO'^ eV shows a small angle clustering within 
the angular resolution of the detectors. The AGASA de- 
tector has found 3 doublets and 1 triplet among 47 de- 
tected events |8]. This analysis is confirmed also by the 
combined data of different detectors |@1 in which 8 dou- 
blets and 2 triplets are found in 92 collected events. This 
evidence can be well understood in terms of a rectilinear 
propagation of protons at the highest energies (E > lO'^ 
eV) with a random arrival of two (three) particles from 
the same source and a source number density of about 
Hs ~ 10-^ Mpc-3 j^. However, the small angle cluster- 
ing may survive in the case of UHE protons propagation 
in IMF 0,111. 

Another remarkable evidence of an almost rectilinear 
propagation of UHECR, in the energy range 2 — 8 x lO'^ 
e V, has been found by Tiny akov and Tkachev 1 1 1 ] . These 
authors have found a correlation between arrival direc- 
tions of UHE particles in the AGASA and Yakutsk de- 
tectors and the directions of several BL-Lac objects, (i.e. 
AGNs with jet directed toward us). The combined evi- 
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FIGURE 1. [Left Panel] Maximum generation energy E™^" defined as min[£g(£,/o), £; 
acceleration energy and Iq is the age of the universe (see text). [Right Panel] Maximal distance rmax 
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where Ei^.^^ is the maximal 



as function of the observed energy E. Tliree merging curves in the left-low corner give r^^^{E) and the dash-dotted curve gives 
rjjj^(£), which numerically is very close to the energy-attenuation length l^tt{E) ■ 
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dences of small angle clustering and correlation with BL- 
Lacs favor a scenario with a quasi-rectilinear propagation 
of protons at energies larger than 10^'^ eV. 



UHECR DIFFUSIVE PROPAGATION 

In order to describe the diffusive propagation in IMF of 
UHECR protons we will use the Syrovatsky 1 12| solution 
to the diffusive equation. Following 0], we will also 
assume a distribution of sources on a lattice; under this 
hypothesis the diffuse flux can be calculated as the sum 
over the fluxes from the discrete sources at distances r,: 



dE, 



exp 



(1) 



4X(E.E„ 
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where b{E) = dE/dt is the proton energy losses (we 
have used b{E) as computed in ifisil ^ and Qi„j — 
{LpK{yg)/El^-„){Eg/EminY^^ is the particle generation 
rate per unit energy with Lp the source luminosity (here 
we assume identical sources with the same luminosity) 
and K{'Yg) ^jg — l (if jg > 2) a normalization constant. 
The function A (/ij/ij,) is 



X{E,Eg) 



de 



m 
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with D{E) the diffusion coefficient. The quantity in 
(|3 describes the squared distance traversed by a pro- 
ton in the direction of the observer, while its energy 



decreases from Eg to E. The Syrovatsky parameter 
X{E,Eg) poses a natural cut to the flux contributing 
sources, it is clear from equation Q that a source at dis- 
tance r > 2y/X iE,Eg) gives a negligible (exponentially 
suppressed) contribution to the flux. 

The Syrovatsky solution ([0 formally includes all 
propagation times up to f — > oo and the generation energy 
of particle is limited only by the maximum energy that 
the source can provide E^nax- Nevertheless, the propaga- 
tion time should be smaller than the age of the universe 
t(), this poses an additional limit to the generation energy 
represented by Ei,{E,tQ), which can be computed evolv- 
ing backward in time the proton energy from E at t = 
up to Eg at t — ffl. Therefore, to limit the propagation 
time, the upper limit of integration £'™^''(£') in equation 
Q is fixed at the minimum between /Smax and Eg{E,to). 
It is interesting to note that at energies E < E^ = I x 10^^, 
where only adiabatic energy losses are relevant, the up- 
per limit of integration in Q is fixed by E^'^^{E) = 
Eg{E,tQ) while at higher energies, where pair-production 
and photo-pion-production energy losses become rele- 
vant, the upper limit of integration is is™^" = iSmax- This 
behavior of £'™''(£') is responsible for the low energy 
steepening of the UHECR diffusive spectrum. In Figure 
n (right panel) we have reported the /s^"" as function of 
the observed energy E, with £max = 

1 X 10^2 eV. 

Let us now concentrate on the diffusion coefficient 
that enters in X{E,Eg) function, see equation Fol- 
lowing J2l, we assume diffusion in a random magnetic 
field with a strength Bq on the coherence length Z^. This 
assumption determines the diffusion coefficient D{E) at 
the highest energies when the Larmor radius of protons 
rL{E) = l(Bo/nG)"'(£'/10'^eV) Mpc becomes larger 
than Ic, namely D{E) = Do{rL{E)/lc)^ with Dq = clc/3. 
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FIGURE 2. [Left Panel] Energy spectrum in the case of So = 1 nG, /e = 1 Mpc and for the diffusion regimes: Kolmogorov 
(continuous line), Bohm (dashed line) and D{E) <x (dotted line). The separation between sources is d = 30 Mpc and the injection 
spectrum index is Yg = 2.7 (see text). The AGASA-Akeno data with the universal spectrum (dash-dotted line) are also reported. 
[Right Panel] The galactic iron-nuclei spectrum computed by subtracting the extragalactic proton spectrum from the Akeno- 
AGASA data. The extragalactic proton spectrum is taken for the case Bq = I nG, !c = I Mpc, d = 30 Mpc, = 2.7 with the Bohm 
diffusion atE < Er. 



At low energy (i.e. when rL{E) < Ic) we have consid- 
ered three different cases: (i) the Kolmogorov diffu- 
sion Dk{E) = DQ{rL{E)/kYI^; (ii) the Bohm diffusion 
Db{E) = DoirtiE) Ik); (iii) an arbitrary case D{E) ^ E" 
with a = 2 as the extreme possibility, also in this case 
D{E) is normalized by Dq at — Ic, that corresponds to 
an energy ~ 10^^ {Bo/nG){lc /Mpc) eV. The diffusion 
length can be evaluated through the interpolation formula 
Idiff = ^d + d/lc, with: A,, = ri and Aj = IdrL/lcY^^ 
in the Bohm and Kolmogorov cases respectively. At dis- 
tances r < ljiff{E) the proton propagation becomes rec- 
tilinear and the fluxes from individual sources on the lat- 
tice are computed in the rectilinear propagation limit. In 
this limit, taking into account the cosmological evolution 
of the Universe, the diffuse flux in our lattice model will 
be: 

" ^ ' ^ (47tE ■ )2 L( £ . ) X (3) 

1 dEg{E,Zi) 
""r^il+zi) dE 

where z,- is the red shift that, according to the standard 
cosmology, is associated to the source at distance r,-, 
the two quantities Eg{E,Zi) and dEg{E,Zi)/dE can be 
computed according to 1 13]. 

Following we fix a reasonable strength of the IMF: 
namely Z?o = 1 nG on the coherence scale = 1 Mpc. 
This choice of (Bq-'c) is compatible with the limits that 
follows from the Faraday rotation measures and with 
a rectilinear propagation regime at the highest energies 
E > lO'^ eV. In order to reproduce, in our lattice model. 



the density of sources as follows from small angle clus- 
tering, we have chosen a separation between sources, i.e. 
a lattice spacing, of c/ = 30 Mpc that corresponds to a 
source space density of = l/d^ — 3.7 x 10^^ Mpc^^. 
In the computation of the fluxes we have assumed an in- 
jection spectrum with a single power law ~ 2.7, start- 
ing from the minimum energy Emm — 1 GeV. 

As follows from equation Q the particles that, as 
a whole, contribute to the diffusive flux are those par- 
ticles produced inside a sphere of radius rmax{E) = 

2^ X{E ,E^^'^), the contribution to the flux is exponen- 
tially suppressed for particles produced at higher dis- 
tances. In Figure ^ (left panel) we have reported, fix- 
ing Ic ~ I Mpc and Bq = I nG the behavior of the 
maximum contributing distance rniax(£^) in the three 
cases Kolmogorov, Bohm, and D{E) °^ E^. The dot- 
dashed line represents the proton attenuation length 
ktt =E{dE/dl)-\ that can be interpreted as the maxi- 
mum contributing distance in the rectilinear propagation 
regime. Figuren(right panel) clearly illustrates the steep- 
ening of the diffusive spectrum at low energies. While 
the energy-attenuation length /att diminishes with energy 
and has the GZK steepening at energy £ ~ 5 x lO'^ eV, 
the diffusive maximum distance rmax{E) increases with 
energy and has a sharp jump at energy isgq — 2 x 10^^ 
eV. The position of this jump is related only to the pro- 
ton energy losses, it does not depend on the magnetic 
field configuration (i.e. on the diffusion parameters). The 
energy Egq corresponds to the proton energy at which 
pair-production energy losses become equal to adiabatic 
energy losses. At this energy the upper limit of integra- 
tion in the Syrovatsky solution changes abruptly from 
Eg{E,t()) to the maximum energy fimax that the source 



can provide. 

Let us consider now the UHECR spectra shown in Fig- 
ure |2l (left panel). In this figure we report the diffusive 
and universal spectra, the latter corresponds to the rec- 
tilinear propagation (see equation 0) in the case of a 
homogeneous distribution of sources (see 1 14| for a de- 
tailed discussion). The effect of the low energy steep- 
ening is clearly seen in the spectra, with the steepen- 
ing energy Es being independent of the diffusive regime 
chosen. According to the results presented in Figure ^ 
(right panel) for ry„ax{E) the flux below the steepening 
energy is largest for the Kolmogorov diffusion and 
lowest in the case D{E) °^ E^, with the Bohm diffusion 
between them. The source luminosity Lp needed to pro- 
vide the observed spectrum is very high, for a distance 
between sources of li = 30 Mpc the required luminosity 
is Lp = 3.0 X 10^^ erg/s. To reduce the required emis- 
sivity one can assume that the acceleration mechanism 
works only starting from a somewhat higher minimum 
energy ^min (see for a detailed discussion). 

It is also worthwhile to stress that the feature of the 
dip, that signals in the experimental data for a proton 
dominated spectrum, is not washed out by the IMF in the 
diffusive approximation. The validity of this approxima- 
tion is related to the validity of the Syrovatsky solution at 
energies E < lO'^ eV. The diffusive equation, and hence 
its solution, are valid only in the case in which the energy 
losses and diffusion coefficient are time independent. At 
energies lower than lO'^ eV, this is not the case, because 
the propagation time of protons approaches the age of the 
universe and the effect of the CMB temperature variation 
with time (red-shift) becomes important. However, our 
computations show a good agreement between the quasi- 
rectilinear regime of propagation and the exact rectilin- 
ear propagation, this agreement shows the approximate 
validity of the Syrovatsky solution at the discussed ener- 
gies. 

Following (3 1^ fisll we shall conclude by discussing 
shortly the transition from galactic to extragalactic cos- 
mic rays. The remarkable feature of the diffusive spec- 
tra is the low-energy steepening at the fixed energy E^ ^ 
1 X lO'^ eV, which provide the transition from extra- 
galactic to galactic CR. This energy coincides approxi- 
mately with the position of the 2nd knee £'2K and gives 
a non-trivial explanation of its value as E2K ~ E^. Like 
in the above-mentioned works we shall assume that at 
£ > 1 X 10^^ eV the galactic spectrum is dominated by 
iron nuclei and calculate their flux by subtracting the 
calculated flux of extragalactic protons from all-particle 
Akeno spectrum. For these calculations we shall fix the 
spectrum with magnetic configuration ( 1 nG, 1 Mpc), the 
Bohm diffusion and a separation between sources on the 
lattice li — 30 Mpc The calculated spectrum of galac- 
tic iron is shown in Figure |2] (right panel) by the dashed 
curve. This prediction should be taken with caution be- 



cause obtained with a model-dependent calculations (as- 
sumption of the Bohm diffusion) and uncertainties in- 
volved in the Syrovatsky solution. However, it is inter- 
esting to note that the iron-nuclei spectrum in Figure 
(right panel) is well described by the Hall diffusion 1 16] 
in the galactic magnetic field at energies above the knee. 

CONCLUSIONS 

We have analyzed the anti-GZK effect in the diffusive 
propagation of UHE protons. This effect consists in a 
sharp increase of the maximum distance rn,ax(£) from 
which UHE protons can arrive (Figure (left panel)). 
The observational consequences of the anti-GZK ef- 
fect is a low-energy steepening of the diffuse spectrum. 
While the shape of the steepening depends on the mag- 
netic field configuration, the steepening energy E^ is 
practically independent of it. The steepening of the spec- 
trum at Es ~ 1 X 10^^ eV coincides with the 2nd knee 
observed in the CR spectra by most of the detectors and 
provides a natural transition from galactic iron nuclei to 
extragalactic protons. 
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